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Abstract 

We show that, unlike other pure b — > d penguin processes, the decay B®(t) — > <pK s 
is dominated by a single amplitude, that of the internal t-quark. The contributions 
of the u- and c-quark operators each vanish due to a cancellation between the (V — 
A) ® (V - A) and (V - A) ® (V + A) matrix elements. Thus, the indirect CP 
asymmetry in this decay probes sin 2/3. Although this cancellation is complete only 
for certain values of the s- and 6-quark masses, the theoretical uncertainty on sin 2(3 
is still less than 10% over most (~ 80%) of the parameter space. By measuring the 
direct CP asymmetry, one can get a better idea of the probable error on sin 2/3. 
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It has been known for many years that the B system is a particularly good place 
to test the standard model (SM) explanation of CP violation. By measuring CP- 
violating rate asymmetries in the decays of neutral B mesons to a variety of final 
states, one can cleanly extract the CP phases a, (3 and 7 This allows one to 
construct the unitarity triangle || and search for the presence of physics beyond the 
SM. 

In the early days of the field, only tree-level decays of B mesons were considered. 
However, it was soon realized that penguin amplitudes could play an important 
role H ffl. For example, the presence of penguin contributions in B d (t) — ► 7r + 7r~ 
can spoil the clean extraction of a (though this can be rectified with the help of an 
isospin analysis jS|). And the clean measurement of 7 via B®(t) — > pK s is completely 
ruined since, for this decay, the penguin amplitude is the dominant contribution. 

Given the importance of such penguin contributions, one is immediately led to 
consider CP violation in pure penguin decays. In the (approximate) Wolfenstein 
parametrization |J of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, there are 
only two matrix elements which have a nonzero weak phase: V td oc exp(— i(3) and 
V u b oc exp(— 27). Thus, assuming that the penguin amplitudes are dominated by an 
internal t-quark, one expects that the b — > s penguin amplitude, which involves the 
product of CKM matrix elements V t bV t *, is real, to a good approximation. Similarly, 
the weak phase of the b — ► d penguin amplitude (VtbV t * d ) is +/3. Knowing that the 
weak phases of B d -B d and B®-B® mixing are, respectively, —f3 and 0, this allows us 
to compute the weak phase probed in various pure-penguin decay asymmetries ||: 

b -> d : Asym(5°(t) -> K°K^) ~ , (1) 

Asym(5°(t) -> (j)K s ) ~ - sin 2(3 , (2) 

b -»• s : Asym(B%(t) -> 0K S ) ~ + sin 2/? , (3) 

Asym(B° s (t) -> 00) ~ . (4) 

The problem with the above analysis is that the 6 — » d penguin amplitude is 
no£ dominated by an internal t-quark. In the quark-level decays b — > uud and 
b — > ccd, the uu and cc quark pairs can rescatter strongly into an ss quark pair, 
giving effective V u bV* d and V c bV* d contributions to the b — ► d penguin decays above. 
Buras and Fleischer have estimated that these contributions can be between 20% 
and 50% of the leading t-quark contribution M. And since the u- and c-quark 
contributions have a different weak phase than that of the t-quark contribution, this 
implies that the weak phase of the b — > d penguin amplitude is not +/3, so that the 
predictions of Eqs. (|l|) and (§) are not valid. On the contrary, due to the presence 
of these several decay amplitudes, one expects that a weak phase cannot be cleanly 
extracted from the measurement of CP asymmetries in pure b —>■ d penguin decays. 
One also expects to observe direct CP violation in such decays. 

Note that these same conclusions do not hold for b —>■ s penguin amplitudes. In 
this case, the CKM matrix-element product associated with the u-quark contribu- 
tion, V u bV* s , is a small fraction (~ 2%) of those of the corresponding c- and t-quark 
contributions, V c bV* s and V t bV t * s . Since the c- and t-quark CKM matrix elements 
are both real, the assumption that the b — > s penguin amplitude is real is a good 
approximation. Thus, the predictions of Eqs. (0) and (ffl) still hold. 
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In this paper, we re-examine the question of the weak phase of the b — > d penguin 
for the exclusive decay B° s {t) — > <pK s [Eq. (0)]. As we will show, although the quark- 
level contributions from it- and c-quarks are non-negligible, at the meson level the 
matrix elements involving the corresponding u- and c-quark operators each vanish, 
to a good approximation, over a large region of parameter space. This is due to 
two factors. First, for a large range of values of the s- and 6-quark masses, it is a 
fortuitous numerical coincidence that the matrix element of the (V — A) £g> (V — A) 
piece of the w/c-quark operator is approximately equal to that of the (V— A)®(V+A) 
piece. Second, because the final state consists of a vector meson and a pseudoscalar 
meson, the full matrix element is proportional to the difference between these two 
pieces. In other words, there is a cancellation between the (V — A) <g> (V — A) and 
(V — A) <g> (V + A) matrix elements. Thus, to the extent that this cancellation is 
complete, the CP- violating rate asymmetry in B®(t) — > <pK s still cleanly probes the 
weak phase (3. 

We begin the discussion by considering the SM effective hamiltonian for hadronic 
B decays of the type b — > dff M: 



W = ^[v fb v; d ( Cl o d lf + c 2 o d 2f ) 



10 

- £ W + V cb V* d cl + v tb v td 4)of] + h.c. (5) 

In the first two terms, / can be a u or a c quark, while in the last three terms, the 
superscript u, c or t indicates the flavour of the internal quark. The operators Of 
are defined as 

Off = d alll Lfpf^Lb a O d 2f = d^Lff-fLb 

Oi 5 = d llM Lbq' llM L(R)q' 0% = d^^Lbpq'^^R)^ (6) 
3 _ 3 _ 

°7,9 = ^d^Lbe^q'^R^q' O° 10 = -d^Lbpe^q'^RiL)^ , 

where R(L) = l±7s, and q' is summed over u, d and s. 0\ and 0% are, respectively, 
the direct and QCD-corrected tree-level operators, 3 „ 6 are the strong gluon-induced 
penguin operators, and operators O7-10 are due to 7 and Z exchange (electroweak 
penguins), and "box" diagrams at loop level. The Wilson coefficients c{ are defined 
at the scale \i ~ m b and have been evaluated to next-to- leading order in QCD. Below 
we give the regularization-scheme- independent values for the c{ for m t = 176 GeV, 
otsijnz) = 0.117, and // = m& = 5 GeV 0: 

Cl = -0.307 , c 2 = 1.147 
4 = 0.017 , c\ = -0.037 , 4 = 0.010 , 4 = -0.045 , 
4 = -1.24 x 10" 5 , c| = 3.77 x 10" 4 , c\ = -0.010 , c\ = 2.06 x 10" 3 , 

4,5 = -<«/#e = P l s/N c , 4 >9 =P e \ 4,10 = 0, l = U,C, (7) 

where iV c is the number of colors. The leading contributions to P % se are given by 
PI = (t) c 2(f + %-M 2 )) and PI = (^){N c c x + c 2 )(f + G{m u fi,q 2 )), where 
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the function G(m, fj,, q ) takes the form 

"7 



/-r/ 2n i /" M m - X(l - X)q 2 

G(m, a,q ) = A / x(l — arjln 5 ax 

Jo zr 



All the above coefficients are obtained up to one-loop order in the electroweak in- 
teractions. The momentum q is the momentum carried by the virtual gluon in the 
penguin diagram. Note that c\ = d 6 , i = u, c, while c\ 7^ c|. This will be important 
in what follows. 

For the decay B — > /1/2, we are really interested in the matrix elements of the 
various operators. We therefore define new coefficients 



-u,c _ ^ r^i / /qn 

c? " </i/ 2 |Oi|S> ' lj 
From Eq. ([/]), the c"' c can be expressed in terms of the function G™'° defined as 
rw.cz \ (/1/2I G(m UjC ,/i,g 2 )Oi (S) 

G - (m *"" ) = (/./,10,1b) • (10) 

One can use models to calculate the functions G"' c , and one finds in general that 
the functions G"' c are process dependent. More importantly, it is a reasonable 
assumption that the functions G™ ,c are independent of i: 

G^ c (m MiC ,/i) = Gf c (m UtC ,fi) . (11) 

This is because the effects of the different Dirac structures of the operators Oj cancel 
in the ratio in Eq. ( |TOD flID| . This implies that the relations between the various c^' c 
are the same as those between the various c*' c (q 2 ). In particular, we have 

4 = 4 ,i = u,c. (12) 

This will be the key ingredient in the analysis below. In the study of non-leptonic 
decays the usual practice is to replace 

-u,c _^ c^ c (q 2 av ) , (13) 
where q 2 v is allowed to vary between m 2 / '4 and m 2 /2 to account for process depen- 
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The structure of the effective Hamiltonian allows us to write the amplitude for 

W s -> <j>K s as 

Af* = A{W S - <j>K s ) = ^(V ub V: d P u + V cb V c * d P c + V tb V t * d P t ) . (14) 

Under the assumption of naive factorization one can write, dropping factors common 
t° Pu,c,ti an d using the fact that 4 = 4 ,i = u,c [Eq. flI2p], 

P u ,c = Ce' C (l - T^) [(Oll) - 2 (O sp )] , (15) 

1 * c 



where 



(Oll) = (^| 57^.(1 " 76)6 B°)(K s \dY(l-l 5 )s\0) , 
(0 5P ) = (0|s(l- 75 )6 B°)<tf fl | d(l + 75)s|0) . 



(16) 

(The operator Osp appears due to a Fierz transformation: (V — A) <g> (V + A) = 
—2(S — P) <S> (S + P).) On the other hand, the contribution from the top penguin 
is more complicated: 

(4 + §-) (0 LL ) + (4 + f ) (Olli) 



+ 



-2(4 + ^) (Osp) + (4 + ^) (Olpi) 
(4 + °f) (Olli) + (4o + ^) (Oll) 



(17) 



where 



(Olli) = (0| s 7 „(l - 7s )s |0) (X s | d 7 "(l - 7s )6 
(0 Lfl i> = (0|^(l + 75 )s|O)(X s | Jy»(l- 75 )6 



(18) 



In the above, we have neglected the contributions from oj$. 
It is convienient to rewrite P UjC and P t as 

= cr(l - 4)* (Oll) , 
P t = a 6 X (0 LL ) + (a 4 - a 6 - ^ai ) (Oll) + (a 3 + a 5 - ^a 9 ) (Olli) , (19) 
where (Olli) = (0 LR1 ), 



Ci + ^, 2 = 4,6,10 
Q + , i = 3, 5, 9 . 



and 



X 



2 (Osp) 
(O 



LL/ 



(20) 



(21) 



It is this latter quantity X which is the focus of our attention in this paper. 
Using the fact that 

(K s (p K )\dl^ ~ 7 5>| 0) = if Ks PKv , 



(22) 



along with the equations of motion for the quarks (we assume that p K = Pd + Ps), it 
is straightforward to show that 



X = 



1-2 



1 



rat 



m b + m s m s + m d 



(23) 
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However, the key point is the following: taking m K = 500 MeV, m;, = 4.9 GeV, 
m s = 100 MeV (all at the 6-quark mass scale), and ~ 0, one finds that X = 0! 
Thus, the matrix elements vanish for u and c but do not vanish for t. The decay 
B° s {t) — ► <pK s is therefore dominated by a single decay amplitude — the t-quark 
penguin contribution — and a measurement of the CP-violating rate asymmetry 
probes the angle j3 [Eq. (Q)]. 

This result is related specifically to two facts: (i) the decay is a pure penguin 
decay, and (ii) the final state <fiK s consists of a vector and a pseudoscalar meson. This 
can be seen as follows. First, consider a decay such as B^(t) — ► 7r + 7r~, which receives 
a tree-level contribution TV u bV* d O T , where the operator O t is of the form (V — A) ® 
(V — A). As we saw above, the matrix element of this operator alone does not vanish. 
Thus, one must consider decays which have no tree-level contributions. That is, only 
processes involving the quark-level transition b — > dss need be considered. Now 
consider a pure-penguin final state containing two pseudoscalars, such as B%(t) — > 
K°K°. Repeating the above analysis, one finds that 



Clearly there is now no cancellation between the two contributions, and one finds 
X ^ 0. Similarly, if the final state consisted of two vector mesons, such as in 
Bd(t) ~* K*°K*°, then the (Osp) matrix element would vanish due to conservation 
of angular momentum, and again we would find 1/0. Thus, only final states 
which consist of a vector and a pseudoscalar can have X = 0. 

In fact, Bg(t) — > 4>K S is the only decay involving a b — > d penguin amplitude for 
which the O u and O c matrix elements vanish. The only other b — ► dss decay whose 
final state consists of a vector and a pseudoscalar is B®(t) — > K°K°*. However, 5° — > 
K°K°* and 5^ — > K°K°* do not factorize in the same way. Specifically, for B® — > 
~K°K *, the is COU pied to the vacuum as in Eq. so that X = 0. However, 
for — > K°K *, it is the which is coupled to the vacuum, which implies that 
the (Osp) matrix element vanishes. Thus, I / for this decay. Therefore it is 
only the decay B®(t) — > <pK s whose indirect CP asymmetry is expected to measure 
sin 2(3. 

Now, the vanishing of the u- and c-quark matrix elements in B®(t) — * <pK s 
depends on two theoretical ingredients: (i) we have assumed that naive factorization 
is valid for this decay, and (ii) we have taken particular values for m s and m&. It 
is important to examine the extent to which these assumptions are justified, and to 
show how the result changes when these assumptions are relaxed. 

First of all, we can see the effect of a nonzero value of X as follows. The 
measurement of the time- dependent rate B®(t) — > <pK s allows one to extract both 
direct and indirect CP-violating asymmetries. These are defined as follows: 



X 



1 + 2 



TTlb — m s Wis + f^d 



(24) 



\A4>K S \2 _ \A<t>K s \2 



I - 



a 



CP 
dir 



\At Ks \ 2 + \Af Ks \ 2 ' 




(25) 



'indir 



\At Ks \ 2 + \Af Ks \ 2 ' 
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where we have assumed that the weak phase of B®-B® mixing is negligible, as is 
the case in the SM. Af Ks is defined in Eq. (|lj), and Af Ks is obtained from Af Ks 
by changing the signs of the weak phases. Using CKM unitarity to eliminate the 
V u bV* d term in Eq. (|HD, Af Ks can be written as 

A*** = ^(V cu e^ + V tu e i5 ^) , (26) 

where we have explicitly separated out the strong phases 5 C and 5 t , as well as the 
weak phase (3. The magnitudes of the CKM matrix elements have been absorbed into 
the definitions of V cu and Vt u - Using this expression for Af Ks , the CP asymmetries 
take the form 



n CP — 



indir vl + + WtuVcu cos ^ cos A 1 J 

where A = 5 t — 5 C . From these expressions, we see that a nonzero value of X 
corresponds to a nonzero value of V cu . This in turn leads to a nonzero value of 
the direct CP asymmetry a^[ r , and also affects the clean extraction of sin 2(3 from 
the indirect CP asymmetry. In order to compute the error on the measurement of 
sin 2/5, we will need to estimate the size of the ratio V cu /V tu , as well as the strong 
phase A. 

The main factor which may contribute to X ^ is if the masses m& and m s 
do not take the respective values 4.9 GeV and 100 MeV. In computing the CP 
asymmetries in Eq. fl2"7|), we will allow m;, and m s to take a range of values. In our 
calculation, we use current quark masses, evaluated at the scale /i ~ m^. For the b- 
quark mass, we take 4.35 < nib < 4.95 GeV. As for the current strange-quark mass, 
there is a great deal of uncertainty in the value of m s . One can obtain m s {n = nib) 
by using as input m s (/i = 1 GeV) and then using QCD running to increase the scale 
up to = nib. However, m s (fi = 1 GeV) is not very well known and can range from 



VPcuVtu sin (3 sin A 
Vl + V 2 CU + 1VtuV cu cos (3 cos A ' 
Vj sin 2(3 + 2V cu Tt u sin (3 cos A 
Vl + VI + 2V tu V cu cos (3 cos A 



0.150 to 0.200 GeV [13]. In addition, it is not clear that perturbative calculations 
are reliable near fi ~ 1 GeV. Given all these uncertainties, we vary m s (nib) in the 
range 0.08 < m s < 0.12 GeV. 

There are also nonfactorizable effects which might give rise to X ^ 0. There 
have been several attempts to calculate corrections to the naive factorization as- 



sumption. One promising approach is QCD-improved factorization [14|, in which 



one systematically calculates corrections to naive factorization in an expansion in 
a s {nib) ~ 0.2 and AQco/^b- Naive factorization appears as the leading-order term 
in this expansion. If we consider QCD corrections to this term, we note that the 
P UjC arise already at 0(a s ), and so they receive no corrections at this order. In 
fact, the P UjC are part of the 0(a s ) corrections to the naive factorization results. 
There are additional a s corrections which can be taken into account by the replace- 
ment (1 + Ti) in Eq. (|19"D , where ~ 0(a s ) are pro cess- dependent 
corrections to the naive factorization assumption. 

Corrections to the P UjC in Eq. ([T5|) may then arise if Eq. (ITT]) does not hold. One 
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can then write 

Pu,c = (1 - ^) K' c (O^) - 2# e (Osp)] = cl\l ~jp)[X + Xi) (Oll) , (28) 



where 



X 1 = • (29) 

c 6 



Note that X\ is complex, so that if one writes X% = \X\\e , one obtains 



|P UjC | = |4' c |(l-— )^/X2 + X 1 2 + 2\X\\X x \cose{0 LL ) . (30) 

In what follows we assume that the effect of X\ is small enough to be absorbed in 
the uncertainty in the quark masses in the expression for X in Eq. (|2"3"D. In other 
words the effect of X\ is essentially taken into account by varying the quark masses 
inX. 

In order to test the robustness of the claim that the indirect CP asymmetry in 
B d (t) — > <pK s measures sin 2f3, we perform the following analysis. We scan the entire 
parameter space, calculating the CP asymmetries a d [ r and af^ dir of Eq. ( f27|) at each 
point in this space. We are especially interested in the quantity 5, which measures 
the fractional difference between the indirect CP asymmetry and the true value of 
sin 2(3: 

r _ a indir ~ 2/3 / Q1 x 

6 = sin 2(3 ■ (31) 

In particular, we wish to compute what fraction of the parameter space leads to a 
given value for 5. This will give us some sense of the extent to which the asymmetry 
in B®(t) — > (f>K s truly probes (3. 

In addition to the masses mj and m s , the CP asymmetries depend on other 
quantities. First, they depend on the momentum transfer q^ v , which we vary between 
m|/4 and ml/2. Second, they involve the CKM parameters p and 77, whose ranges 
are taken from Ref. [|HJ. Third, one also needs values for the d- and c-quark masses. 
We take m d (fi = m b ) = 6 MeV and 1.1 < m c < 1.4 GeV @. Finally, note that, 
while P UjC depend only on the matrix element (Oll), Pt a ls° depends on (Olli) 
[Eq. (P^DI- Using factorization, one can write 



(O ll1 ) UF? S ( 



ml' 



(O ll ) f Ks At{m? K ) ' 
where the various formfactors and decay constants above are defined by 
{K s {q)\di^-lh)u\ 0) = if Ks q» , 
(0(9, e) \s lfll5 b\ Sj(p)) = (Mz + M^Af 



(32) 



e*-(p-g ), 
(P - <?) 2 



^ e*-(p~g) 
2 M B + ikL 



M 2 - Ml 
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(0 \&fs\ e)) 



(p ~ if 



(K s (q) \uj^b\ B° d (p)) = F{ 



K S 



(P + <?f " 



B 



m 



k s 



(p - qY 



+F K ° 



(p - q) 2 



(p - qY 



Ma, 



(33) 



Using = 0.237 GeV and the various form factors calculated using the light 
cone sum rules [1171 we obtain r' w 1.5. 



We are now in a position to calculate the CP asymmetries and CL^ dir in 
Bg(t) — > 4>K S [Eq. (|27|)1) as well as the deviation of the indirect CP asymmetry 
from the true value of sin 2(3. In Fig. [I], we plot the fraction of the parameter space 
for which \8\, the error on sin 2(3 as extracted from the indirect CP asymmetry 
[Eq. fl3ll) 1, is less than a given value, |5| maa; . From this figure, we see that sin 2(3 
can be obtained with an error less than 30% over virtually the entire parameter 
space. And this error is reduced to about 10% in 80% of the parameter space. 
While this should not be interpreted statistically as some sort of confidence level, 
it does indicate that it is quite likely that (3 can be extracted from the indirect CP 
asymmetry in B®(t) — > 4>K S with a rather small error. 
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Figure 1: The fraction of the parameter space for which the error \5\ in extracting 
sin 2(3 from af^ dir is less than a given value, |5| mox . 
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Of course, if the time-dependent rate for B®(t) — ► <pK s is measured, we will have 
more information than just af^ dir : we will also measure the direct CP asymmetry 
a d [ r . Since a d [ r vanishes if X = 0, its value may help us determine the extent to 
which affair really measures sin 2f3. At first glance, the correlation between a d ( r and 
a indir appears airtight: if a d [ r is found to vanish, then this must imply that X = 0, 
so that a^dir yields sin 2f3. Unfortunately, things are not quite so straightforward: 
a d [ r is also proportional to the strong phase difference A [Eq. (|2"T|)]. Therefore, if 
A ~ 0, then a d ? will vanish even if X ^ 0. Thus, this possibility must be taken into 
account in evaluating the correlation between the measurements of af^ dir and a d f r . 




0.025 0.05 0.075 0.1 0.125 0.15 0.175 0.2 



I dir max 

Figure 2: The fraction of the parameter space for which a d [ r is less than a given 
value, |o^| ma:r . 

In Fig. |2|, we plot the fraction of the parameter space for which a d £. is less than 
a given maximal value. As is clear from the figure, we expect a d [ r to be at most 
20% - larger values would point to the presence of new physics. Fig. |3| shows, for 
a given value of |a^>|, the fraction of the parameter space for which \5\ is less than 
a particular maximum error (|#| maa; = 5%, 10%, 15%, 20%). From this plot, we see 
that if a d [ r is measured to be 0.1, one can obtain sin 2/5 from af^ dir with an error of 
5% (20%) in ~ 55% (~ 95%) of the parameter space. If a d [ r is found to be tiny, 
then this is probably due to the fact that X ~ 0, since 5 < 5% over ~ 90% of the 
parameter space. However, as discussed above, this does not hold over the entire 
space since a d [ r can be small if A ~ 0, while 1^0. 

There is one technical point which is worth mentioning here. Naively, one would 
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Figure 3: The fraction of the parameter space for which 5 can be measured with a 
maximum error of \S\ max = 0.05, 0.1, 0.15 and 0.20, as a function of \a^\. 



expect that, as a%. increases, the error 5 will also increase. That is, one expects 
that the fraction of the parameter space with a given maximum error |#| maa; should 
decrease with increasing a^[ r . This is because the prediction that the indirect CP 
asymmetry in B° s {t) — > 4>K S measures sin 2(3 depends on the fact that X = 0, and 
a c J* r increases with increasing X. However, Fig. |3] does not behave in exactly this 
way: although the fraction of the parameter space with a given \S\ max does indeed 
roughly decrease with increasing a^, this decrease is not monotonic. For example, 
the curve for \S\ max = 0.05 flattens out at around \a%?\ ~ 0.04, and even turns up 
before falling again for ^0.1. The explanation for this behaviour is as follows. 
From Eq. (|27|), we see that \a^[ r \ depends, among other things, on the strong phase 
difference A and the ratio r = V cu /Vt u - As q\ v is varied, keeping other parameters 
fixed, we find that, as r increases, so does sin A. However, \5\ depends on cos A, 
which obviously decreases as sin A increases. Thus, there is a region of parameter 
space where, as r and sin A increase, the error \S\ actually decreases due to the 
presence of the cos A term. It is this effect which is responsible for the flattening 
out and the slight rise in the curves in Fig. |3|. 

Finally, one can take a most conservative point of view, and ask what is the 
maximum error |<5| maa: over the entire parameter space for a given measurement of 
l a cfa>l- This is shown in Fig. [|. Depending on the value of \a^f r |, \5\ max is between 
30 and 40%. As was the case in Fig. |, the downturn of the curve as |a^| increases 
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Figure 4: The maximum error |<5| max for the entire parameter space as a function of 
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is due to the dependence of 5 on cos A. 

Of course, if X does indeed vanish, this may have some negative practical im- 
plications. Specifically, since there are fewer contributions to the amplitude for 
Bg — > 4>K S , one might suspect that the branching ratio will be smaller than that of 
other pure b — » d penguin decays. This is indeed what we find: taking p = 0.18 and 
rj = 0.36 |15||, P u = P c =0, and using the form factors in Ref fl~T|l , we obtain 



BR[B° S -> <f)K s 



BR[B° d 



P+7T 



0.003 . 



(34) 



P+7T ] 



30 x 10" 6 pi we find BR[B° -> (j)K s 



Using the measured BR[B® 
10 -7 , which is very small. 

Fortunately, the above analysis for B®(t) — > <f)K s also applies to the decay 
Bg(t) — > <p(1680)K s , where 0(1680) is a radially excited 0. We expect the branch- 
ing ratio for B° s — > 0(168O)iT s to be almost a factor 10 larger than B° s — > <pK s |]19 . 
This is because the form factor for B® — > (or, in general, for any i? — > light me- 
son) probes the high-momentum tail of the wavefunction. As the radially excited 
0(1680) has more high-momentum components, the form factor for B° s —>■ 0(1680) 
is enhanced relative to B® — > 0. 

Note that the measurement of CP violation in B®(t) — > (j)K s probes (3 [Eq. (fj)]. 
If the measurement of j3 as extracted in this mode disagrees with the measurement of 
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(3 from Bj(t) — > J/^K S , this will indicate the presence of new physics in the b — > s 
penguin amplitude, i.e. in the b — > s flavour-changing neutral current (FCNC) pOf . 
Similarly, the value of (3 extracted in B° s {t) — ► <pK s or B®(i) — > (b(1680)K s can 
be compared with that found in B®(t) — > J/^K S . Assuming that the weak phase 
of B®-B® mixing is tiny — and this can be tested by looking for CP violation in 
Bg(t) — ► J/^4>, for example — a discrepancy between these two values points clearly 
to new physics in the b — > d FCNC. This new physics might affect B®-B® mixing 
and/or the b — > d penguin amplitude. Now, it is quite likely that CP violation in 
Bg decays can only be measured at hadron colliders, since one needs an extremely 
large boost in order to resolve the rapid B®-B® oscillations. Since hadron colliders 
produce copious amounts of B% and B® mesons, it should be possible to perform the 
B®(t) — > 4>K S and B®(t) — > 4>K S analyses simultaneously, since the final state is the 
same. Thus, by measuring (3 in these decay modes, at hadron colliders one can test 
for the presence of new physics in both the b — ► s and b — > d FCNC's. 

To summarize: in general, b — ► d penguin decays receive contributions from 
internal u, c and t quarks. Because of this, one cannot cleanly extract information 
about weak phases from measurements of CP-violating rate asymmetries. In this 
paper, we have shown that this does not hold for the decay B®(t) — > 4>K S . Due to 
a fortuitous cancellation between the (V — A) <g> (V — A) and (V — A) <g> (V + A) 
matrix elements, the contributions from the u- and c-quark operators each vanish. 
Thus, Bg(t) — > 4>K S is dominated by a single decay amplitude, that of the internal 
t-quark, so that the indirect CP asymmetry measures sin 2/3. This is the only b — ► d 
penguin decay for which this occurs. 

Of course, this cancellation is not rigourous: it depends on our choosing partic- 
ular values for the s- and 6-quark masses. However, we have shown that, over most 
(~ 80%) of the theoretical parameter space, the difference between the indirect CP 
asymmetry and the true value of sin 1(3 is less than 10%. Furthermore, by using 
information about the direct CP asymmetry, one can get a better handle on the 
probable error on sin 2/3. For example, should a^[ r be found to be tiny, this increases 
the likelihood that sin 2(3 is being extracted with a small error: for a^f r ~ 0, the 
error on sin 2(3 is at most 5% over almost 90% of the parameter space. 

Because the u- and c-quark contributions to B®(t) — ► <pK s are very small, one can 
expect that the branching ratio for this decay will be reduced. This is indeed what 
is found: we estimate that BR[B® — > (j)K s \ ~ 10~ 7 . However, the same analysis also 
applies to the decay B° s (t) -> 0(1680)7^, where 0(1680) is a radially excited 0. The 
branching ratio for this decay is expected to be about 10 times larger than that for 
B° s (t) ^ <f>K 3 . 

Finally, we note that B®(t) — > <pK s will probably be studied at a hadron collider. 
These same experiments can simultaneously measure the CP asymmetry in B^(t) — > 
<f)K s , which also probes sin 2(3 in the SM. By comparing the CP asymmetries in these 
two modes with that measured in B®(t) — > J /ipK Sl one can look for the presence of 
new physics in both b — > d and b — ► s transitions. 
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